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The fermionic and bosonic sectors of the 2-site Hubbard model have been exactly solved by 
means of the equation of motion and Green's function formalism. The exact solution of the t-J 
model has been also reported to investigate the low-energy dynamics. We have successfully searched 
for the exact eigenoperators, and the corresponding eigenenergies, having in mind the possibility 
to use them as an operatorial basis on the lattice. Many local, single-particle, thermodynamical 
and response properties have been studied as functions of the external parameters and compared 
between the two models and with some numerical and exact results. It has been shown that the 
2-site Hubbard model already contains the most relevant energy scales of the Hubbard model: the 
local Coulomb interaction U and the spin-exchange one a consequence of this, for 

some relevant properties (kinetic energy, double occupancy, energy, specific heat and entropy) and 
as regards the metal-insulator transition issue, it has resulted possible to almost exactly mime the 
behavior of larger systems, sometimes using a higher temperature to get a comparable level spacing. 
The 2-site models have been also used as toy models to test the efficiency of the Green's function 
formalism for composite operators. The capability to reproduce the exact solutions, obtained by 
the exact diagonalization technique, gives a firm ground to the approximate treatments based on 
this formalism. 



PACS numbers: 71.10.- 



71.10.Fd 



o 

> 
O 

o 



I 

o 
o 



X 
J3 



I. INTRODUCTION 

Two are the aspects that gave so much popularity to 
the Hubbard model: the richness of its dynamics that is 
thought to permit a description of many puzzling issues 
like metal-insulator transition, itinerant magnetism, elec- 
tronic superconductivity, and the simplicity of the Hamil- 
tonian structure that let one speculate about the possi- 
bility of finding the exact and complete solution for any 
realization of the underlying lattice. Anyway, although 
the model has been studied more than any other one in 
the last fifty years, very few exact results are available 
and what we have mainly regards either finite clusters or 
the infinite chain (i.e., the ID case). For finite clusters 
of 2^'^ or 42-3.4,5,6 gi^eg \i jg possible to find the complete 
set of eigenstates and eigenvalues of the Hamiltonian and 
compute any quantity by means of the thermal averages. 
However, it is not easy at all, although possible in prin- 
ciple, to extract valuable and scalable (i.e., which can be 
used to find the solutions of bigger and bigger clusters 
and, ultimately, of the infinite lattice cases) information 
regarding the effective excitations present in the system, 
the operators describing them and their dynamics. For 
the infinite chain neither, we have all the information we 
wish; the Bethe Ansatz is a very powerful tool, but is 
severely limited as regards the range of applicability of 
the self-consistent equations it supplies and the quanti- 
ties for which it gives an answer. 

In this manuscript, in order to overcome the limita- 
tions discussed above, we have exactly solved the Hub- 
bard model, on a 2-site cluster, completely within the 
equations of motion and the Green's function formalism. 
By using this approach, we have had the possibility to 
find the complete set of eigenoperators of the Hamilto- 



nian and the corresponding eigenenergies. This informa- 
tion has been really fundamental as it permitted a deeper 
comprehension of the features shown by the properties 
we have analyzed. It is worth noting that the Hubbard 
model on a 2-site cluster is the smallest system where 
both terms of the Hamiltonian (i.e., kinetic and electro- 
static) are effective and contributes to the dynamics. 

By properly tuning the value of the temperature, we 
have found that the 2-site system can almost perfectly 
mime, as regards relevant properties such as the kinetic 
energy, the double occupancy, the energy, the specific 
heat, the entropy and fundamental issues such as the 
metal insulator transition, the behavior of larger clusters 
and of the infinite chain. The tuning of the tempera- 
ture is necessary in order to get a comparable effective 
level spacing (bigger the cluster, lower the spacing), i.e. 
to excite the correct levels: the relevant energy scales 
are present although the relative positions of the lev- 
els are affected by the size of the system (only two k 
points!). The very positive comparisons with exact re- 
sults (Bethe Ansatz, exact diagonalization) and numer- 
ical data (quantum Monte Carlo, Lanczos) support the 
idea that a lot of physics can be described and under- 
stood within this very small system, for which there is 
the possibility to know the analytic expressions for all the 
quantities under study. As regards the relevant scales of 
energy, this 2-site system has demonstrated to contain 
all the necessary ingredients to describe many features 
coming from the strong electronic correlations and also 
appearing in the lattice case. Two of the three relevant 
energy scales, which are thought to be present in the 
Hubbard model, naturally emerge: the local Coulomb 
interaction U and the spin-exchange one J = ^ that 
is, in principle, extraneous to the original purely electro- 
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static Hamiltonian and is dynamically generated by the 
combined actions of the two terms of the Hamiltonian. 
As useful guide to better understand the low-energy dy- 
namics we have also solved the t-J model and presented 
the solution in parallel with the one found for the Hub- 
bard model. 

This analysis, which have resulted to be really relevant 
by itself as we have got a much better understanding of 
some energy scales and internal parameter dynamics, has 
been worth to be performed also as a prelude of the lattice 
analysis. In fact, the eigenoperators we have found, both 
in the fermionic and the bosonic sectors, can be used in 
the lattice case as a basis for the Green's functions. In the 
strongly correlated systems, the interactions can alter so 
radically the dynamics of the original particles that these 
latter lose completely their own identities^. Actually, 
some new objects are generated by the interactions and 
dictate the physical response. They are not so easy to be 
identified: their number, exact expression and relevance 
can only be suggested by the experience and, when avail- 
able, by exact and/or numerical results. For instance, 
one can choose: the higher order fields emerging from 
the equations of motion, the eigenoperators of some rele- 
vant interacting terms, the eigenoperators of the problem 
reduced to a small cluster, etc. In the last years, we have 
been focusing our activity on the study of strongly corre- 
lated electronic models like Kondo, t-J , p-d, Hubbard by 
means of the Composite Operator MethodS^^iiSiii'^ that is 
based on two main ideas: one is the use of composite fields 
as basis for our Green's functions, in accordance to what 
has been discussed above, and the other one is the ex- 
ploitation of algebra constraints (e.g., the Pauli principle, 
the particle-hole symmetry, the Ward-Takahashi identi- 
ties, ...) to fix the correct representation of the Green's 
functions and to recover the links among the spin and 
charge configurations dictated by the symmetries. It is 
worth noticing that the Composite Operator Method is 
exact in itself. An additional approximation treatment 
is needed when we deal with large or infinite degree-of- 
freedom systems; in this case we have to treat in an ap- 
proximate way the otherwise intractable hierarchy of the 
equations of motion generated by the projection proce- 
dure. If no approximation is necessary (finite and reason- 
ably small degree-of- freedom systems) , the COM cannot 
do else than give the exact solution. According to this, 
the COM gives the exact solution also for the two systems 
under analysis in this manuscript: the two-site Hubbard 
and t-J models. Whenever, instead, we should resort 
to an approximate treatment to close the hierarchy of 
the equations of motion generated by the projection pro- 
cedure, we expect some limitations connected with the 
chosen approximation. For instance, we get only the first 
moments correct if we truncate the equations of motion 
hierarchjiiS. On the other hand, it is really worth not- 
ing that we properly take into account: the interaction 
term of the Hamiltonian by using as basis operators its 
eigenoperators"'^^; the short-range correlations by using as 
basic fields the eigenoperators of the problem reduced to 



a small cluster^'^; the presence of a Kondo- like singlet at 
low-energy by properly closing the equation of motion of 
an ad hoc chosen composite operator—. 

Obviously, we are aware that the exact diagonalization 
of this very small system takes less than one afternoon 
to any graduate student. Then, the reader could wonder 
why we decided to study such a system so in detail. Well, 
the reasons are many and some have been already pointed 
out above: 

• Any graduate student can surely compute eigen- 
states and eigenvalues of these Hamiltonians in one 
afternoon, but the analytic computation of Green's 
functions and correlation functions, in terms of the 
former eigenstujj, is not that straightforward as one 
can think. At the end of the day, the time saved 
in computing eigenstates and eigenvalues instead 
of eigenoperators is almost fully recovered if you 
also take into account the time needed by the com- 
putation of the physical properties (Green's func- 
tions and correlations functions). Then, the pos- 
sibility to have scalable information putting, alto- 
gether, almost the same effort become quite tempt- 
ing for anyone. At any rate and for the sake of com- 
pleteness, we report in Appendix the expressions 
of Green's functions and correlations functions in 
terms of eigenstates and eigenvalues. 

• The knowledge of the exact eigenoperators of a sys- 
tem is invaluable as they could be used as correct 
starting point for the application of the projection 
methods^^iiSiiLiSiiS to strongly correlated systems 
whose minimal model^"'^ is the exactly solved one. 

• The Green's function formalism for composite oper- 
ators is extremely complicated-^. The 2-site Hub- 
bard model can be used as toy model to fully 
explore this formalism and evidence the difficul- 
ties connected with the treatment of composite 
operators with non-canonical commutation rela- 
tions. Within this system, the appearance of zero- 
frequency functions can be safely handled and re- 
solved. The links among the different channels 
(fermion, charge, spin, pair) can be studied in de- 
tail. 

• The 2-site Hubbard model, according to its sta- 
tus of minimal modeUi, contains the main scales 
of energy related with the interactions present in 
the Hamiltonian. The exact solution in terms of 
eigenoperators permits to individuate which are the 
composite fields responsible for the relevant transi- 
tions. These latter can be then used to efficiently 
study larger clusters. 

• The absence of the three-site term, for obvious ge- 
ometrical reasons, in the derivation of the 2-site 
t-J model from the 2-site Hubbard one permits to 
push further the study of the relations between the 
two models. It is possible to exactly individuate 
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the low-energy contributions and study them sepa- 
rately. 

• Last, but not least, the capability of miming the nu- 
merical results for larger clusters (sometimes tun- 
ing the temperature and, consequently, the effec- 
tive level spacing) opens the possibility to provide 
a testing-ground for the numerical techniques. 

In the following, we define the models, give the self- 
consistent solutions in the fermionic and bosonic sectors 
and study the local, single-particle, thermodynamic and 
response properties of the systems; the eigenoperators, 
eigenstates and eigenvalues of the system are also given 
and analyzed in detail. 



II. THE 2-SITE HUBBARD AND t-J MODELS 

The Hamiltonian of the Hubbard model^" for a A^-site 
chain reads as 

H^Y. - ^^5^J)c\{)c{]) + (2.1) 



/I is the chemical potential, c^{i) = ^c|(i), c|(i)^ is the 

electronic creation operator at the site i in the spinorial 
notation, U is the on-site Coulomb interaction strength, 
n„{i?) = cJ.(i)co-(i) is the charge density operator for spin 
a at the site i and 



tij — —2taij — —t ^ cos [k {i — j)] a{k) 



(2.2) 



where t is the hopping integral, a is the lattice constant, 
aij is the projection operator on the nearest-neighbor 
sites and a{k) — cos(fca). In the momentum representa- 
tion, the kinetic term Ht of the Hamiltonian 1)2. 1() reads 
as 



Ht = -2t^Q(fc)ct(fc)c(fc) 



where k assumes the values with / = 0, 

Na ' 

for periodic boundary conditions and 



(2.3) 

,iV- 1 

, N for open boundary conditions. We will study a 
2-site cluster within periodic boundary conditions: N = 
2 and fc = 0, ^. 

The Hamiltonian of the t-J modelSiiSSiS^ for the same 
cluster and boundary conditions reads as 

H = Y,[-2ta{k)~^i]^\k)ak) 

k 

where ^{i) = [1 — n(i)] c(i) is the fermionic composite 
operator describing the transitions n = ^ n = 1, 
= is the total charge {fi — 0) and spin 



(/i = 1, 2, 3) density operator at the site i, = (1,1?), 
cr^ — (—1,(7) and (7 are the PauH matrices. 

We will extensively use the following definition for any 
field operator ^'(i) 



vE'"(z) = Ea,,vI/(j) 



(2.5) 



In particular, for the 2-site system we have 5*" (0) = 
* (a), (a) ^ (0) and *"'(i) = (*"(i))" = *(i). 



III. THE FERMIONIC SECTOR 
A. The equations of motion and the basis 

1. The Hubbard model 

After the Hubbard Hamiltonian [Eq. (|2.1|) ]. the elec- 
tronic field c{i) satisfies the following equation of motion 



i^c(^) =-Mc(^)-2te"(^) + C/77(^) 



(3.1) 



with rj{i) ~ n{i)c{i) — —^ai.nk{i)c{i). According to this, 
we can decompose c{i) as 



(3.2) 



where ^(i) and 77(1) are the Hubbard operators and de- 
scribe the transitions n — ^ n — 1 and n = 1 n = 2, 
respectively. Moreover, they are the local eigenoperators 
of the local term of the Hubbard Hamiltonian and de- 
scribe the original electronic field dressed by the on-site 
charge and spin excitations. They satisfy the following 
equations of motion 



= - 2te"(z) - 2t7r(z) 

i^7?(z) = iU-fiMi) + 2tnii) 



(3.3) 



with 



7^(^) = 7rCT''rv(0c"(«) +e(»)c^"(«)r/(*) 



(3.4) 



where 71.^(1) = c^(i)(T^c(i) is the total charge [fi = 0) and 
spin (/i = 1, 2, 3) density operator at the site i in the 
Hubbard model. We use a different symbol to distinguish 
it from the analogous operator, z/^, we have defined in the 
t-J model. The field Tr{i) contains the nearest neighbor 
charge, spin and pair excitations dressing the electronic 
field c{i). 

The field 7r(i) can be also decomposed as 



(3.5) 
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where 



(3.6) 



These latter, which are non-local eigenoperators of the 
local term of the complete Hubbard Hamiltonian, satisfy 
the following equations of motion 



i^77,(z) = iU- ^i)7Js{i) + 2t7^{i) + 2tC{i) 



(3.7) 



By choosing these four operators as components of the 
basic field 



( ^(0 \ 

6(0 



(3.8) 



we obtain a closed set of equations of motion. In the 
momentum space, we have 



i^V(fc) = e(fc)VXfc) 
where the energy matrix e(fc) is 



e(fc) = 



( -[i~2ta{k) ~2ta{k) 
U-^i 
At 
\ 2t 



(3.9) 



-2t -2t \ 

2t 2t 
-fi + 2ta{k) Ata{k) 
2ta{k) U-fi J 
(3.10) 



It is worth noting that ip(i) is an eigenoperator of 
the Coulomb term of the Hubbard Hamiltonian for any 
lattice structure according to the local nature of the 
interactioni^. 



2. The t-J model 

For the t-J model we can obtain a closed set of equa- 
tions of motion by choosing as basic field 



where 




(3.11) 



(3.12) 



We have 

i^V(fc) = e(fc)V'(fc) 



with 



s{k) = 



-^-2ta{k) -2t + 2Ja{k) 
-^ + 2ia(fc)-4J 



(3.13) 



(3.14) 



It should be noted that the t-J model exactly repro- 
duces the Hubbard model in the regime of very strong 
coupling (i.e., U ^ t). The three-site term, appearing in 
the derivation of the t-J model from the Hubbard one, 
is completely absent in the 2-site system for obvious ge- 
ometric reasons and the given Hamiltonians are exactly 
equivalent in the strong coupling limit. In particular, 
we note the following limits: lim[/^t^s(i) = C(i) and 
lim;7>t r]{i) = lim;7>f r]s{i) = 0. 



B. The Green's function 

1. The Hubbard model 

Let us now compute the thermal retarded Green's func- 
tion G{k,uj) ^ T{TZ U)]) that, after the equa- 
tion of motion of ijj{k) (|3.9|) . satisfies the following equa- 
tion 



[w-e(fc)]G(fc,w) =/(fc) 



(3.15) 



where I{k) = ^|^(i), '0t(j)}^^ ^ is the normalization 
matrix. The subscript E.T. means that the anticommu- 
tator {. . .} is evaluated at equal times. and Tl[. . .] 
stand for the Fourier transformation and the usual re- 
tarded operator, respectively. (. . .) indicates the thermal 
average in the grand canonical ensemble. The solution of 
Eq. 13.151 by taking into account the retarded boundary 
conditions, is the following 



G{k,uj) 



4 



r(") 



(fc) 



uj - En{k) -f i(5 



(3.16) 



The spectral weig htS cr(")(fc) can be computed by 
means of the following expression 



(3.17) 



where A is a matrix whose columns are the eigenvectors 
of the energy matrix e(fc). 

The energy spectra i?„(fc) are the eigenvalues of the 
energy matrix e(fc) 



Ei{k) = -n-2ta{k) 
E2{k) ^ -fj.-2ta{k) + U 
Esik) = + 2to(fc) - 4Ju 
Ei[k) =-11 + 2ta{k) +AJu + U 



(3.18) 
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+ 64t2 — U). It is worth noting 
^ which is the value of J we get 



where Ju = ^ {VU 

that limu^t Ju - -[j 
from the derivation of the t-J model from the Hubbard 
one in the strong coupling regime. One can easily check 
that the energy spectra En{k) represent the energy of the 
single-particle transitions among the eigenstates given in 
App.EI 

The normalization matrix I{k) has the structure 



m 



hi 





/l3(fc) 








I 22 







hz{k) 





/33(fc) 






V 



(3.19) 



/24(fc) /44(fc) / 

The explicit form of its entries are 
hi 



h2 — 



n 
2 



h3{k) = A + a{k){p-l22) 
hA{k) = -A-a(k)p 
l33{k) = -2{p-l22)-2a{k)A 
hi{k) = h2 



(3.20) 



with 



n — {nil)) 

A = (^(^)et(^))-(r?"(^)r/t(^)) 



(3.21a) 
(3.21b) 

(3.21c) 
(3.21d) 
(3.21e) 



A gives a measure of the difference in mobility be- 
tween the two Hubbard subbands. p contains charge, 
spin and pair correlation functions. In order to com- 
pute the electronic Green's function Gcc — [cc'''] ) — 
Gil + 2012 -I- G22 we only need the following quantities 



a«(fc) = /n 




t{k)Ii3{k) 



^^ihk) 



af^k) 



1 



a(fc)/i3(fc) 1 



(3.22) 



,(3) 



(fc) 
(fc) 



0^^{k) = 

-a^^ik) 



-/l3(fc) 



2t 



(3.23) 



-^22 

^(2) 
-^22 



(fc)=0 



(fc) 



-^22 
1 

2' 



-a(fc)/i3(fc) 



42^(fc) = -^«(fc)^13(fc) 



4Jr, 



(3.24) 



.(4)^1.^ _ 



1 



(fc)---a(fc)/i3(fc)(l- 



U + 8Ju 



Equation (|3.1t)|) does not define uniquely the Green's 
function, but only its functional dependenceii. The 
knowledge of the spectral weights cr*^"^ (fc) and the spec- 
tra En{k) requires the determination of the parameters 
A and p, together with that of the chemical potential 
/i. The connection between the parameter A and the el- 
ements of the Green's function and the equation fixing 
the filling allow us to determine these two parameters as 
functions of the parameter p 



2(1 - Gil - 

' ^11 ~ '-'22 



C22) 



(3.25) 



where the correlation functions Cab = \jpa{'i)^l{'i)^ and 

Cab = ^V"" (*)V'&(*)^ can be computed, after the spectral 
theorem, by means of the following equations 

"^^ (3.26) 



k n—1 



with 



T„{k) = tanh 



2T 



(3.27) 



The presence of parameters related to bosonic corre- 
lation functions within the fermionic dynamics (e.g., the 
parameter p) is characteristic of strongly correlated sys- 
tems. In these systems, the elementary fermionic exci- 
tations are described by composite operators whose non- 
canonical anticommutation relations contain bosonic op- 
erators. Then, according to the general projection pro- 
cedure, which is approximate for an infinite (i.e., with 
infinite degrees of freedom) system and coincident with 
the exact diagonalization for a finite system, the energy 
and normalization matrices generally contain correlation 
functions of these bosonic operators. 

The parameter p could be computed through its def- 
inition (see Eq. I3.21c|l . In this case we need to open 
the bosonic sector (i.e., the charge, spin and pair chan- 
nels), to determine the eigenoperators and the normal- 
ization and energy matrices and to require the complete 
self-consistency between the fermionic and the bosonic 
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sectors. Moreover, we run into the problem of fixing the 
value of the zero-frequency constants—. This procedure 
will be widely discussed in the next section where the 
bosonic sector and the relative channels will be opened 
and completely solved. However, we can use another pro- 
cedure to fix the parameter p without resorting to the 
bosonic sector. We have to think over the reason why 
the parameters fj,, A and p appear into our equations: we 
have not fixed yet the representation of the Hilbert space 
where the Green's functions are realized. The proper 
representation is the one where all the relations among 
operators coming from the algebra (e.g., — ng-) and 
the symmetries (e.g., particle-hole, spin rotation, ...) are 
verified as relations among matrix elements. In the case 
of usual electronic operators the representation is fixed 
by simply determining the value of the chemical poten- 
tial. In the case of composite fields, which do not satisfy 
canonical anticommutation relations, fixing the represen- 
tation is more involved and the presence of internal pa- 
rameters (e.g., fi, A and p) is essential to the process 
of determining the proper representation. The require- 
ment that the algebra is satisfied also at macroscopic 
level generates a set of self-consistent equations, the local 
algebra constraints^ that will fix the value of the internal 
parameters^. 

In particular, for the 2-site system the equation 



(3.28) 



coming from the algebra constraint f(i)?7^(j) — 0, to- 
gether with Eqs. 13.251 constitute a complete set of self- 
consistent equations which exactly solves the fermionic 
dynamics allowing to compute the internal parameters 
(i.e., fi, A and p) for any value of the model (t, U) and 
thermodynamical (n, T) parameters. This procedure is 
clearly extremely simpler than that requiring the open- 
ing of the bosonic sector. It is worth mentioning that 
the system of self-consistent equations can be analyti- 
cally solved as regards the parameters A and p as func- 
tions of the chemical potential /i. It is possible to show 
that the equation for this latter parameter, as function of 
the internal (i.e., model and thermodynamical) parame- 
ters, exactly agrees with that coming from the thermal 
averages (cfr. Eq. IA42|I as it should be since the model 
is exactly solved. This further confirms the validity, the 
effectiveness and the power of the method used to fix the 
representation. 

In the next section, we will see that also in the case 
of the bosonic sector, the determination of the proper 
representation will be obtained by means of the local al- 
gebra constraints which will fix the value of the zero- 
frequency constants. It will be also shown that the two 
procedures for computing the parameter p are equivalent, 
as it should be, and give exactly the same results for the 
fermionic dynamics although with remarkably different 
effort. It is worth noticing that the use of the local alge- 
bra constraints, which is unavoidable to fix the represen- 
tation both in the fermionic and in the bosonic sectorsii, 
permits to close the fermionic sector on itself without re- 



sorting to the bosonic one also in the lattice caseS^ where 
a fully self-consistent solution of both sectors, although 
approximate, is very difficult to obtain. 



2. The t-J model 

In the t-J model, the energy spectra (i.e. 
values of £(fc)) read as 

Ei{k) = -fi-2ta{k) 
E2(k) = -^i + 2ta{k) - 4:J 



the eigen- 



(3.29) 



In particular, E2{k) corresponds to the single-particle 
transitions between the single occupied states and the 
2-site singlet. The latter is the ground state at half fill- 
ing. C(fc) drives the corresponding fermionic excitation 
between the single occupied states and the 2-site singlet 
and the appearance of the most relevant scale of energy 
at low temperatures J. 

The normalization matrix I{k) has the following en- 
tries 



hi 



1 



n 
2 



h2{k)^C^,+a{k) (^^X"- 2 

/22(fc) - -2/i2(fc) 



(3.30) 



In the t-J model the spin and charge correlator x" is 
defined as = S^=o ('^^i agreement with 

the strong coupling nature of the model. The spectral 
weights cr'-"-' (fc) have the following expressions 



a[\\k) 
a^^^ik) 



-lh2ik) 



'12 
.(2) 



'22 

,(2) 
'22 



{k) = 
(fc) = 





/l2(fc) 



(fc)-o 

(fc) -/22(fc) 



(3.31) 



(3.32) 



(3.33) 



Also in this case, as for the Hubbard model, we could 
compute x" opening the bosonic sector for the t-J model. 
Again, following the reasoning given in the previous sec- 
tion, a simpler and completely equivalent procedure relies 
on the exploitation of the algebra constraint C(*)C^(*) — 
—2({i)£^^i). According to this, the parameters fi and x" 
can be computed by means of the following set of self- 
consistent equations 



77 = 1 
C22 — 



-2Ci2 



(3.34) 
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IV. THE BOSONIC SECTOR 
A. The equations of motion and the basis 

1. The Hubbard model 

a. The spin and charge sectors. After the Hubbard 
Hamiltonian [Eq. (|2.1() ]. the charge (/i = 0) and spin 

= 1, 2, 3) density operator n^{i) — c\i)a^c{i) satisfies 
a closed set of equations of motion, which describes the 
spin and charge dynamics in the system under study, once 
we choose as basic field 



where 

























- wlii) 















(4.1) 



where 
/o(*) 



:ct(j)a^c"W-c"t(j)a^c(z) 

---vH^)v{^)-4{i)d^{i) 



(4.2) 
(4.3) 
(4.4) 
(4.5) 

(4.6) 
(4.7) 



The components of suggested by the hierarchy of 

the equations of motion, are either hermitian or anti- 
hermitian as densities or currents should be at any order 
in time differentiation. (j>^{i) satisfies the following equa- 
tion of motion in momentum space 

d 



LUb{k) 





^4i[l -a(A:)] 







-2t \ 

[/ 

U 2t 

U 2t 

8i 

8t y 



(4.9) 



b. The pair channel. Within the analysis of the dy- 
namics of the Hubbard model, another relevant bosonic 
operator is the pair operator p{i) = C|(z)c|(i). The set 
of composite fields 



P{^) = 



( P(i)(z 

p(3)(j 

p(5)(, 
y p(6)(,; 



(4.10) 



where 

p(i: 

p(2 
p(3 
p(4 
p(5 
p(6. 



p{i) 

ct(z)cj(z) 

ctW<W + ^t(*)c?(*) 

p(i)n°'{i) 
p{z)v''Hi)r{i) 



(4.11) 
(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 

satisfies the following closed set of equations of motion 
. d 



(4.17) 



(4.8) where 



-2t [l + a{k)] 







-2t [l + a(fc)] 
-2^1 








U 

U-2fi 
-2t[l + a(fc)] 





2. The t-J model 

The basis for the bosonic sector in the t-J model is 
jiven by 







-2t [l + a{k)] 
U-2fi 









u 



2([/-/i) 
-2t [l + a{k)] 



-2t[l + a(fc)] 
U-2^ j 



r 



for the charge channel, and by 



Vk{i) 



(4.18) 



(4.20) 
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for the spin channel. We have defined 



where the relevant entries of the normalization matrix 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



In the momentum space we have the following closed 
sets of equations of motion 



d 
d 

^Ol'fikiq) = LJs{q)(pk{q) 



(4.25) 
(4.26) 



with 



^^ciq) = 



-2t 



-4t[l-a{q)] 



(4.27) 



and 



i^siq) 



/ -2t -2j\ 

-U[l~a{q)] -At[l-a{q)] 
2J 

\ 4J[l-a(g)] / 

(4.28) 

B. The Green's function 



InRef.IH we have shown that the retarded and causal 
Green's functions contain substantially different informa- 
tion according to the unavoidable presence of the zero 
frequency constants (ZFC). In particular, we have re- 
ported on the relations between different types of Green's 
functions and on the correct order in which they should 
be computed. According to this, in the bosonic sector, 
we have to start from the causal Green's function and 
not from the retarded one, as we correctly did in the 
fermionic sector. 



1. The Hubbard model 

a. The spin and charge sectors. After Eq. 14.81 
the thermal causal Green's function G^{k,u}) = 
T (T [0/j(«)0ji(j)] ) satisfies the following equation 



/,,(fc)=.F([0^(»),</.t(j)]) 



E.T. 



are 





= 




(4.30) 


/l2p(fc) 


= 4[1- 




(4.31) 


/l3M(fc) 


= 




(4.32) 




= 8[1- 




(4.33) 




= -[1- 


- a(fc)]q(^) 


(4.34) 


/l6M(fc) 







(4.35) 



with 



C^, = (c"(z)ct(z)) 
g(°) = 16d 

3 

fc=i 



(4.36) 
(4.37) 

(4.38) 
(4.39) 



T[. . .] stands for the usual time ordering operator. 
The solution of Eg. 14.291 is the following oneii 



G^(fc,w) = ~2iTTT^{k)S{uj) 

6 

E' C 

i=l 



uj — Eif,{k) + iS uj — Eii,{k) — i6 

(4.40) 



where 

XT [0{U - tj) ((/.^(z)<^t (j)) - 0{t, - U) (0t (J>^(^))] 

(4.41) 

is the zero frequency function and is undetermined at 
this level unless to compute the quite anomalous Green's 
function appearing in its definition, which involves anti- 
commutators of bosonic operators^l. The general defini- 
tion of a zero frequency function in terms of eigenvectors 
and eigenvalues is given in Eq. IA12I The primed sum is 
restricted to values of i for which Eu,{k) ^ 0. 

The Eib{k) are the eigenvalues of the energy matrix 
UJb (fc) 



Euik) 

E2b{k) 

EM 
Eib{k) 

EM 
Eebik) 



-2iV2[r 
2tV2[l - 

-U - AJu 

-AJu 

AJu 

U + AJu 



(4.42) 

(4.43) 
(4.44) 
(4.45) 
(4.46) 
(4.47) 



[uj - LOb (fc)] (fc, u) = I^{k) 



(4.29) 



According to this, the primed sum in Eg. I4.4()l does not 
contain the i — 1 and i = 2 elements for fc = (i?i,26(0) = 



9 



0) and the zero frequency function r^(fc) reduces to the a and b = 1, 2, computed through Eq. 13.171 have the 
constant = r^(0). The spectral weights ^[^(fc) with following expressions 



,(2) 



(k) 



^l-a(fc) U {tq^''^ [1 + a{k)] + 4[7 [1 - a{k)] C^^} + 2C" js [1 + a{k)f - [1 - a(/fc)]|] 



V2 {S [1 + a{k)f - t/2 [1 - a(/c)]} 
4i!7 [1 - a{k)] [4{U + Uu)C^2 - tq^"^] 



{U + 8Ju) {16<2 [1 + a{k)] + 2U{U + 4Ju)} 

AtU [1 - a{k)] [16 JyCfa + t<?^^'] 
" ([/ + 8J[/) {16i2 [1 + a{k)] - WJu} 

I 



(4.48) 
(4.49) 

(4.50) 



and 



(fc) 



(fc)=- 



Ejbjk) (i) 



(4.51) 



2t 



(fc) 



^zb(fc) 

2i 



'11/^ 



(fc) 



We note the sum rules 

('■) 



6 

E 



(4.52) 



(4.53) 



In order to finally compute the Green's function 
Giifj, {k,uj) we should fix the internal parameters x", d 
and the zero frequency constant Fn^. 

The parameter is directly connected to the Green's 
function. Let us consider the correlation function 
(fc,w) = ^ {4'fj.{'i)4'li{j)) which is linked to the causal 
Green's function through the spectral theorem 



(fc, c^) = - (l + tanh ^) 3 [G^ (fc, lo)] 



27rr^(5(cj) + 27r 



J2 S[cj-E,,ik)]- 

i=l 



(fc) 



Then, 



X? = 3(0«(z)4^)"(z) 



g-/3£;,6(fc) 
(4.54) 



(4.55) 



The computation of the parameter d requires instead 
the opening of the pair sector. Otherwise, it could be 
computed through the following relation with the param- 
eter p 



1 



P 



(4.56) 



already given in the fermionic sector (see Ea. l3.21c|l . once 
we note that x" = X? + (4^H«)'/'o^^"(0^ 



The zero frequency constant Fn^ cannot be directly 
connected to any correlation function at this level and 
its determination requires the use of local algebra con- 
straints. In particular, we can use the following relation 



(4.57) 



where the double occupancy D — (ri|(j)n|(i)) is given by 

D = I22 — C22- 

Eauations l4 . 5 51 H . 561 and [4 . 5 71 constitute a complete set 
of self-consistent equations which allow to compute the 
parameters x" > d and Fn^ and then to determine the 
Green's function Gn^ {k,uj). 

It is worth noting that, once C"^, Cfj, P and D are 
computed in the fermionic sector, the bosonic correla- 
tion functions can be easily obtained. We have the fol- 
lowing expressions for the relevant correlators and zero- 
frequency constants 



.n-2D-2a 
1 + 26 



(4.58) 



X? = ('/'|,'^(^)4'^"(*))=« + 2i?-2a 



1 



3(71 -2D- 2a) 
1 + 2& 

3(n -2D- 2a) 



Fii 



4(1 + 6) _ 
1 



l + 2b 



n + 2D - 2a - Ap 

(4.59) 

n + 2D - 2a - Ap 

(4.60) 



1 



[2n + AD-2a 



+b{n + 2D-3 
F113 



n-2D-2a 



1 + 26 
2(n-27:>)(l + 6)-2a 
1 + 26 



Ap 



(4.61) 
(4.62) 
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where 



where 



a = -(C,",-2Ci"2)coth 



2t 
T 



b=- 



U + 8Ju 



coth 



U + 8Ju 



■ coth 



2T 



coth 



2T 



2T 

U + Uu 
U + 8Ju 



coth 



(4.63) 

2T 
(4.64) 



Taking into account the following local algebra con- 
straints 

D,2o{h^)^{n{t)9H^)) = ^C:, (4.65) 
^220 ih «) = (.9(*).9t(*)) - 2n - - 4d 

= 2{n -Ad- 2p) = 2{n - x" + 2p) (4.66) 

Di2^ {t, i) = {n^{i)9l{i)) = 2C^, (4.67) 
^223 (*, *) = (53(*)gl(*)) - 2n + ix" - Xc + 4d 



2(n 



2p) 



(4.68) 



rp(fc) = lim w 

X T [e{t - t') {P{r)P\j)) - 9{t' - t) (pt(J)P(^))] 

(4.80) 

The primed sum is again restricted to values of i for which 

Erp{k) ^ 0. 

The Eip{k) are the eigenvalues of the energy matrix 

LOp (k) 



Eip{k) ^ -2^1 + U - 2t[l + a{k)] 
E2p{k) = -2/i + [/ + 2t[l + a{k)] 

E^p{k)^-2^l+]^[U -Q{k)] 
Eip{k) ^ -2^l+]^[W - Q{k)] 
E^p{k)^-2^l+]^[U + Q{k)] 
EQp{k) = -2n+'^[W + Q{k)] 



(4.81) 
(4.82) 

(4.83) 
(4.84) 
(4.85) 
(4.86) 



and performing similar calculations we can also obtain 
(we omit the expressions of T22^l for the sake of brevity) 



(n"(z)gt(^)) = _2C,", 
{g^{i)g\i))^-2{n-Ad-2p) 
Tl2^, = 



(4.69) 
(4.70) 
(4.71) 



h. The pair channel. After Eq. 14.171 the thermal 
causal Green's function Gp{k,uo) = T {T [P{i)P\j)\) 
satisfies the following equation 



[uj - LUp (k)] Gp (fc, Lu) = Ip{k) 



(4.72) 



where the relevant entries of the normalization matrix 
Ip{k) = T{[P(i),P\])\)^^^^ are 



hip{k) 


= 1 — n 


(4.73) 


Ii2p{k) 


= [l + a(k)]C'^, 


(4.74) 


hspik) 


= Ii2p{k) 


(4.75) 


hipik) 


^n-Xc + "^daik) 


(4.76) 


hbpik) 


^2[l + a{k)]C'^2 


(4.77) 


Imp{k) 


= 2D -2-1 + 2da{k) 


(4.78) 



where 7 — {n°'{i)ni(i)ni{i)) . 

The solution of Eq. 14.721 is the following oneii 



Gp{k,uj) = -2mTp{k)S{uj) 



4'\k) 



-f3LJ 



1 



-/3a; 



i-^ — Eip{k) +i5 LU — Eipik) — \5 



(4.79) 



withQ(fc) = + 16t2[l -ha(A;)]2. The£;,p(fc) are zero 
only in isolated points of the parameter space (n, T, U) 
(see Eas. 14.811 - 14.8611 . According to this, the zero fre- 
quency function Tp{k) is identically zero except in these 
points, where it could be finite. In this treatment, for the 
sake of simplicity, we neglect these isolated points. 

The spectral weights <Jiip{k), computed through 
Eq. 13.171 have the following expressions 



<^iipi^) = ^[^i3p(fc) + hAp{k) - Ii5p{k) - hepik)] 



(4.87) 



'^iipik) = ^hh3p{k) + lupik) + /i5p(fc) - hepik)] 

(4.88) 



+ ^^HrJ^miipik) - 2/i4p(fc) + /i6p(fc)] 



4Q(fc) 



Q{k) + U 



AQ{k) 



(4.89) 
(4.90) 



Q{k) 
Q{k) + U 
4Q(fc) 



2/np(fc)-2/i4p(fc)+/i6p(fc)] 



(4.91) 



(6) t[l + a(fc)] ^ Q{k)-U ^ 

(4.92) 
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In order to finally compute the Green's function 
Giip{k,LL!) we should fix the internal parameter 7. 
The determination of the parameter 7 requires the 
computation of another Green's function GDik,Lu) — 
T {T[n{i)D{j)]) , where D{i) — ri|(j)n|(i). Otherwise, 
we could resort to the following local algebra constraint 



The general solutions of Eqs. 14.941 are the following 



Gc {q,(^) = -2'nrrc{q)5{uj' 
2 U) 



{p{i)p^{i)) = 1 - n + D 



(4.93) 



i=l 



1 



It is worth noting that within the pair channel we can 
also compute the parameter d directly from its definition: 
d — {p"' {i)p^ [i)) . This finally opens the possibihty to 
compute fully self-consistently the fermionic and bosonic 
sectors at once. It can be shown that the two procedures 
(the fully self-consistent one and the one presented at 
length in the previous sections) give exactly the same 
results, as it should be according to the exact nature of 
the proposed treatment. 



2. The t-J model 

After Eqs. 14.251 the thermal causal Green's func- 
tions Gc{q,uj) = T {T [ip{i)Lp'^ {j)]) and Gsk{q,uj) = 



Gsk {q, t^) 

4 M 



uj~Etc{q) + ^S Lu ~ E.ciq) - 

(4.104) 



i=l 



uj - (q) +iS LU - Eis (q) - iS 

(4.105) 



where 



^c{q) = 7; 1™ ^ 

Z uj—>-0 



X :f [e{t t') (^(*)^t(j-)) „ e^^t' - 1) (v^^jX^)] 

(4.106) 

r.(<z) ^ 



lim Lo 

2 w^o 



X T 



satisfy the following equations 



e{t ~ t') [M^^lU)) - e{t' - 1) (^I(j>fe(* 

(4.107) 



[uj - LUc (q)] Gc [q, uj) = Ic{q) 
[uj - ujs (q)] Gsk {q, = Isiq) 



(4.94) 
(4.95) 



The primed sum is again restricted to values of i for which 
E^cAk) + 0. 

The Eic (q) and Eig (q) are the eigenvalues of the energy 
matrices udq) and tOsiq): respectively 



where the relevant entries of the normalization ma- 
trices Iciq) = J" {[^{i),tpHj)])j^,^^ and /^(g) = 



^{ Vk{i),^i{j) are 




E2c{q) 
and 


^ -2tV2[l 




/lie (9) = 


(4.96) 




= 2tV2[l - 




/i2c(g)-4[l-a(g)]Cri 


(4.97) 


Eis{q) 




/22c (g) = 


(4.98) 


E2s{q) 


= -2iV2[l 


- 


and 




Essiq) 


= 2JV2[1 - 





= 

= 4[1 
= 



a{q)] C^, 



hisiq) 
h2s{q) 
hzs{q) 

hisiq) 

with 



fe=l 



(4.99) 
(4.100) 
(4.101) 

(4.102) 



(4.103) 



Eisiq) 



(4.108) 
(4.109) 

(4.110) 
(4.111) 
(4.112) 
(4.113) 



According to this, the primed sums in Eas. l4.lOH do not 
contain the elements for q = and the zero frequency 
functions Tc^siq) reduce to the constants Tc = rc(0) and 
Tg = rs(0), respectively. The spectral weights acl{k), 
computed through Eq. 13.171 have the following expres- 
sions 



^^2'c(9)='^i2i(9)-2Cri[l-a(g)] 



.(2) (^^ - 



The spin rotational invariance makes Gsk independent 
from the index k. According to this and for the sake of 
simplicity, we have omitted in the text the index k in the 
expressions of the related quantities: Is, Es, as and F^. 



42c (9) 
411(9) 
4?l(9) 



-'y^liq) 



/l- 


a{q) 


{q)\ 




[1- 


a{q)Y 


/i- 


a{q) 



(4.114) 
(4.115) 
(4.116) 
(4.117) 
(4.118) 
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In order to finally compute the Green's function 
Gc{q,oj) we should fix the zero frequency constant r^. 
This latter cannot be directly connected to any correla- 
tion function at this level and its determination requires 
the use of local algebra constraints. Let us consider the 
correlation function Dc{q,uj) = J- (^ip{i)(p^ {j)') which is 
linked to the causal Green's function through the spec- 
tral theorem 



D,{q,uj) = -{1 + tanh ^) 5 [Gc {q, w) 



2.r,5(c.) + 2nY! S[u.~ E^q)] ^„,^,^(,) 

(4.119) 



Then, we can use the following relations in order to com- 
pute Tc 



We have the following expressions for the relevant zero 
frequency constants and correlators 



2t 

Tn, = n + coth -C^, 

Tl2c = 
T22C = 

Tn, = 3n + 3 coth ^C^^ + coth 
X? = n + 2coth|cfi 



(^"(z)nt(z)) ^ -2Cfi 

Of 

(n"(z)nt(z)) = 4coth-Cfi 

^ _ 3n + 6coth§$Cfi _ 



3x? 



1-2 coth ^ 1-2 coth ^■ 



(4.127) 

(4.128) 
(4.129) 

(4.130) 

(4.131) 
(4.132) 
(4.133) 

(4.134) 



Diic = {v{i)v{i)) = n 

D22C (*,*) = (n(z)nt(*))^2(n-Xc) 



(4.120) 
(4.121) 
(4.122) 



Xc is directly connected to the Green's function 

X? = (4.123) 

In order to finally compute the Green's function 
Gils (fc,w) we should fix the internal parameter x" and 
the zero frequency constant Fn^. 

The parameter x" is directly connected to the Green's 
function. Let us consider the correlation function 
Dsk {q,^) ~ ^ (^Vk{'i)^\{j)^ which is linked to the causal 
Green's function through the spectral theorem 



UJ \ 



Dsk ((?, w) = - ( 1 + tanh — j 3 [Gsk (g, w)] 



2^r,(5(w) + 2ttJ2'H^- ((?)] - 



ai'\q) 



Then, 

X? = 3(j.3(*)i^3"W> 



(4.124) 



(4.125) 



The zero frequency constant Ts cannot be directly con- 
nected to any correlation function at this level and its de- 
termination requires the use of local algebra constraints. 
In particular, we can use the following relation 



Diisk {hi) = {vs.(i)i''i[i)) = n 



(4.126) 



Equations 14.1251 and 14.1261 constitute a complete set 
of self-consistent equations which allow to compute the 
Green's function Gii^fe (q,w). 



V. THE RESULTS 

In the previous two sections, we have given a detailed 
summary of the analytical calculations which lead to 
computable expressions for internal parameters, zero fre- 
quency constants, correlation and response functions for 
the fermionic and bosonic (charge, spin and pair) sec- 
tors of the Hubbard and t-J models. In this section, we 
analyze the related results for the local, single-particle, 
thermodynamical and response properties. Hereafter, t 
will be used as energy scale. 



A. The A and p parameters 

The A and p parameters are computed by solving the 
system of self-consistent Equations 13.251 and 13.281 The 
behavior of A, as function of the filling n, is shown in 
Fig- Q (top panel). A is defined as the difference be- 
tween the hopping amplitudes computed between empty 
and singly occupied sites (i.e., {£.{i)^H'i))) and between 
singly and doubly occupied sites (i.e., (r/(i)?7^(i))). In 
this sense, A gives a measure of the electron mobility. 
The following features have been observed: the hopping 
amplitude coming from ^ (77) prevails below (above) half- 
filling and is always negative [see inset Fig. ^ (top panel)]; 
the absolute value of A diminishes on increasing the tem- 
perature T and on decreasing the Coulomb repulsion U. 
Below half filling, the electrons move preferably among 
empty and single occupied sites; the same happens to 
the holes above half filling. This explains the prevalence 
of one hopping amplitude per each region of filling. In 
addition, by decreasing the temperature or by increas- 
ing the Coulomb repulsion we can effectively reduce the 
number of doubly occupied sites present in the system 
and, consequently, favor the mobility of the electrons and 
increase the absolute value of A. These are also the rea- 
sons behind the filling dependence: its absolute value 
first increases with the number of available almost free 
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FIG. 1: (top) The A parameter as function of the filhng n 
for T = 1 and (7 = 0, 2, 4, 8, 12 and 16 (inset: i and r] 
components for T = 1 and (7 = 4); (bottom) Cfi as function 
of the fiUing n for T = 2 and J = 0.05, the ED data (4 x 4) 
are for T = and from Ref. I2rI 

moving particles (n < 0.5), then decreases when the num- 
ber of particles is such to allow some double occupancies 
and to force the reduction of the number of empty sites 
(0.5 ^ f^ ^ 1) and finally vanishes when, in average, the 
sites are all singly occupied (n = 1); the behavior above 
half filling is equivalent and is related to the holes. 

The behavior of CJ^ within the t-J model, as function 
of the filling n, is shown in Fig^ (bottom panel). Be- 
low half filling and for rather high values of the Coulomb 
repulsion [/, A and exactly coincide as the t] com- 
ponent of A completely vanishes. According to this, all 
the features discussed above for the parameter A are also 
observed for Cfj^ (J acts like 1/C7). Anyway, we have to 
report a really small dependence on the exchange inter- 
action J that is fully effective only at n = 1 where the 
mobility is already zero as no empty site at all is left. It is 
worth noting that the exact result for the 2-site system at 
T = 2 almost exactly reproduces the Exact Diagonaliza- 
tion (ED)25 data for a 4 x 4 cluster at T = 0. This shows 
that, by increasing the temperature in a 2-site system it 
is possible to mime a cluster of bigger size at a lower tem- 
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FIG. 2: The p parameter as function of the fiUing n for T = 
and T = 2 at (7 = and (7 = 4. 

perature, at least as regards some of its properties. This 
can be understood by thinking to the level spacing in the 
two cases (i.e., bigger the cluster lower the spacing) and 
to the value of temperature needed to excite those levels. 
Clearly, not all the properties are in a so strong relation 
with the relative level positions, but depend on the ab- 
solute positions of them. Anyway, this also shows that 
the relative energy scales/levels are already present in the 
2-site system. We have used a value of J twice smaller 
than the one used within the numerical analysis accord- 
ing to the difference, between the 2-site system and any 
larger system, regarding the value of the exchange energy 
J appearing in the derivation of the t-J model from the 
Hubbard one2&. Actually, the derivation is pathological 
just for the 2-site system and gives a J larger of a factor 
2. 

As regards the -p parameter, in Fig |21 we report its 
dependence on filling n for two values of the on-site 
Coulomb repulsion ?7 = and U — \ and for two tem- 
peratures T = and T — 2. At low temperatures and for 
n < 1.26, the value of the p parameter is mainly negative 
according to the antiferromagnetic nature of the spin cor- 
relations in the ground state; at higher temperatures the 
antiferromagnetic spin correlations get weaker and the 
sign of the parameter changes. The filling n, the temper- 
ature T and the Coulomb interaction U rule the balance 
between the spin, the charge and the pair correlations. 

B. The chemical potential and the density of states 

The chemical potential fi is one of the self-consistent 
parameters and has been computed by solving the system 
of self-consistent Eg nations 13.251 and 13.281 As a function 
of U the chemical potential increases for all studied val- 
ues of temperature T and filling n. In particular, we 
have the maximum increment for U < 20; for higher 
values of U the chemical potential saturates. A similar 
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FIG. 3: The chemical potential jj, as function of the filling n 
(top) for i7 = 4 and T = 0.01, 0.1, 0.2, 0.3, 0.4, 0.5 and 1 and 
(bottom) for f/ = 4, 12 and T = 1, 2. The FTL data (4 x 4) 
are taken from Ref. l27l 
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FIG. 4: The thermodynamic A^t(m) and the traditional N{uj) 
{n — 1 and S = 0.1) densities of states as functions of the 
scaled chemical potential jj, — U/2 and the frequency uj, re- 
spectively, for T = 2 and f/ = 0, 8, 16 and 32. 



behavior has been reported by numerical data on larger 
clustersSSi22422i2£. In Fig. 13 we show the filling depen- 
dence of the chemical potential for some values of the 
temperature T at U = 4 (top panel) and for U = 4, 12 
and T = 1, 2 (bottom panel). The Finite T emp erature 
Lanczos (FTL) data (4x4) are taken from Ref. 1271 At low 
temperatures, the chemical potential has a step-like be- 
havior which recalls the discreteness of the energy levels. 
At higher temperatures, the step-like feature is smeared 
out due to the thermal hybridization of the energy levels. 
It is again really noteworthy the agreement between the 
2-site results and the numerical data. No temperature 
adjustment has shown to be necessary in this case as the 
temperature is already quite high [cfr. Fig. 1121 (bottom)]. 
Generally, the agreement is better for high temperatures 
and high values of the Coulomb repulsion. Anyway, there 
is a competition between the temperature smearing ef- 
fects on the level spacing and the possibility to access the 
states with energy Jjj at lower temperatures for higher 
values of U {Jjj oc 1/C/). Actually, the number of acces- 
sible states is comparable for medium temperatures and 



too large for very high temperatures. This should explain 
the deviations for J7 = 4 and T — 1, and U = 12 and 
T = 2 at low doping. 

By looking at these plots we could be induced to con- 
sider the possibility of a metal-insulator transition driven 
by the Coulomb repulsion U and controlled by the tem- 
perature T. Obviously, no such transition is possible in 
a finite system that is always in a metal paramagnetic 
state for any finite or zero value of the Coulomb repul- 
sion U ; we will show, in Sec. IV HI that the Drude weight 
is always finite. Such tricky behavior of the chemical po- 
tential can be understood by thinking at the nature of the 
grand canonical ensemble we are using for our thermal 
averages. As we are very far from the thermodynamic 
limit, there is no equivalence at all among the micro-, 
the grand- and the canonical ensembles. In particular, 
whenever we speak about temperature and chemical po- 
tential in a finite system we have just to think in terms 
of mixtures of quantum mechanical states with different 
energies and numbers of particles. To get a deeper com- 
prehension of this issue we can define a thermodynamic 
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N{uj) is computed as follows 
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FIG. 5: The thermodynamic Nrif^) and the traditional N{lj) 
{n — 1 and 5 — 0.1) densities of states as functions of the 
scaled chemical potential fj, — U/2 and the frequency ui, re- 
spectively, for (7 = 8 and T = 0.02, 0.4, 0.8 and 2. 



density of states NT{fi) through the following equation 



(5.1) 



which is a generalization to finite temperatures of the 
identical relation existing, at zero temperature, between 
the filling n and the usual density of states N{uj). In gen- 
eral, the two densities of states coincide only at zero tem- 
perature and for systems, like the non-interacting ones, 
whose density of states is independent on the chemical 
potential /i. From the definition H5.1() . we can simply 
compute the thermodynamic density of states A^t(a*) by 
differentiating the filling n with respect to the chemical 
potential fi [Nt{p) = dn/ dfi]. It is worth mentioning that 
the thermodynamic density of states Nt{p-), apart from a 
factor is simply the compressibility of the system. A 
vanishing value of this quantity denotes the impossibility 
for the system to accept more particles and, obviously, 
the presence of a gap. The thermodynamic density of 
states NxifJ-) enormously facilitates the comprehension 
of the chemical potential features that are singled out by 
the differentiation procedure. The usual density of states 



N{lo) 



(fc) 



(5.2) 



n=l k 



Both densities of states are reported in Figs. 01 and [S] as 
function of the chemical potential (/x) and the frequency 
(uj), respectively. We have computed N{uj) by substi- 
tuting, for obvious graphical reasons, the Dirac deltas 
iS{uj)) with Lorentzian functions (^ ^a^^a ) with S = 0.1. 
In reporting the thermodynamic density of states NT{fJ.) 
we have scaled the chemical potential /i by its value at 
half filling {U/2) as we want to make the comparison for 
this particular value of n. While N{oj) always presents a 
gap at a; = (the overlapping tails are due to the finite 
d and the non-zero kinetic energy ensures a metallic be- 
havior), Nt{p-) shows a gap only above a certain value 
of the Coulomb repulsion for a fixed T (see Fig. and 
below a certain temperature for a fixed U (see Fig. O . 
We want to emphasize that, at high temperature, Nt{p) 
is capable to mime the behavior expected for the lattice 
system. In particular, NtIp) presents two well defined 
structures (i.e., the two Hubbard subbands) that continu- 
ously separate on increasing J7 (see Fig. 0}. This behavior 
is the one expected, after many analytical and numeri- 
cal results, for the bulk at dimension greater than one 
and it is a realization of the metal-insulator transition 
according to the Mott- Hubbard mechanism. It is worth 
noting that, for this small system, the critical Coulomb 
repulsion, at which the gap appears, is a function of the 
temperature T. 

The density of states N{u)) is also reported in order to 
give an idea of the positions of the poles and of the rel- 
ative intensity of the spectral weights, although we had 
to cut the highest peaks in drawing the pictures. We can 
observe two relevant features. In the non-interacting case 
{U = 0) only two poles/peaks are present; those coming 
from the band Ei{k) and separated by the bandwidth 
W — At. On increasing the Coulomb potential U the 
three scales of energy present in the system clearly man- 
ifest themselves: the exchange interaction Jjj, the band- 
width W and the Coulomb repulsion U. The eight poles 
group in two main structures separated by U. Within 
any structure the four poles are separated according to 
the combination of ±2t {bandwidth separation) and the 
presence or absence of the exchange interaction Jjj. The 
bandwidth W is obviously U independent and generates 
a rigid two peak structure for any band En{k). On the 
contrary, the exchange energy Ju decreases on increas- 
ing U with a consequent reduction of the resolution of 
the peaks within the Hubbard subbands. The other rel- 
evant feature is the redistribution of the total spectral 
weight on increasing the temperature. At low tempera- 
tures the poles/peaks that do not contain the exchange 
interaction Jjj have negligible spectral weights; the 2-site 
singlet is the ground state at half-filling. On increasing 
temperature the total spectral weight redistributes and 
the other poles get more and more spectral weight due 
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FIG. 6: The double occupancy D: (top) as function of the 
fiUing n for [/ = 1, 2 and 4 at T = 1, the qMC data (12 x 12) 
are for T = 1/6 and from R,ef. l3ll: (bottom) as function of the 
Coulomb repulsion U for n = 0.5, 0.75 and 1 at T = 1, the 
BA data are for T = 0. 



to the thermal hybridization of the energy levels. 



C. The double occupancy D 

In the non-interacting case [U — 0) we have D = 
At zero temperature, the double occupancy D vanishes 
for n < 0.5; below this value of filling we have, in average, 
less than one electron in the system and only for a finite 
temperature we can get some contributions by states with 
a finite double occupancy. At half filling we have the 
following exact formula for the double occupancy at zero 
temperature 



D 



2Ju 



U + 8Ju 



(5.3) 



The filling dependence of the double occupancy D (top 
panel) is reported in Fig. El for several values of the 
Coulomb repulsion U at temperature T = 2. There are 
also reported some quantum Monte Carlo (qMC) data^ 



for a bigger cluster (12 x 12). As in the case of C^i within 
the t-J model, we note that the exact results for the 2- 
site system at temperature T = 1 very well reproduce the 
quantum Monte Carlo (qMC) data for the bigger clus- 
ter at T = 1/6. The same explanation obviously holds. 
A comparison with the exact results from Bethe Ansatz 
(BA) are also reported in Fig. |H| (bottom panel). Also 
in this case the 2-site system manage to reproduce the 
bigger cluster data (actually, the Bethe Ansatz (BA) sys- 
tem is the ID bulk) by increasing the temperature. The 
discrepancy at half filling can be understood as a conse- 
quence of the difference in the definition of the exchange 
energy J discussed in the Section regarding the internal 
energy. Obviously, the discrepancy is larger where the 
exchange interaction is mainly effective (i.e., at half fill- 
ing and for intermediate-strong values of the Coulomb 
interaction) . 



D. Thermodynamics 

1. The energy E, the specific heat C and the entropy S 

The energy per site E is computed as thermal average 
of the Hamiltonian divided by the number of sites. Its 
filling dependence is reported in Fig. [7| (top panel) to- 
gether with some Exact Diagonalization (ED) data for 
bigger clusters^^. As expected, the energy E increases 
as the Coulomb repulsion U increases. For n < 1 and 
high values of the Coulomb repulsion the kinetic term 
prevails on the Coulomb one. The opposite behavior is 
observed for n > 1. The behavior as a function of the 
temperature T and Coulomb repulsion U are reported in 
Fig. El (middle and bottom panels, respectively) together 
with some data coming from numerical analysis for big- 
ger clusterg-^'9^'^^. By using a higher temperature for 
the 2-site data is again possible to get an extremely good 
agreement with the numerical results [see Fig. (top and 
bottom panels)] . In the case of the temperature behav- 
ior instead, the agreement is obtained without tuning any 
parameter [see Fig. [7| (middle panel)]. 

The discrepancy at half filling and low temperatures 
[see Fig. \7\ (middle panel)] is a consequence of the dif- 
ferent value of the exchange energy J appearing in the 
derivation of the 2-site t-J model from the Hubbard one 
(see detailed discussion about the behavior of CJ\ re- 
ported in Fig. 13 (bottom panel). Sec. IV A|) . According 
to this, at half-filling and low temperatures we have also 
reported the 2-site results for U = 24 and J = 1/6 (t- 
J model) and obtained the expected agreement. The 
reason why no effect is evident for lower fillings and high 
temperatures is that the exchange interaction is really ef- 
fective only at half-filling (in average one spin per site is 
necessary) and low temperatures (the fiuctuations should 
be small). 

In Fig.|Slwe show the temperature dependence of the 
specific heat C = dE/dT for several values of U at 
n = 0.75 and n = 1. For n 7^ 1 we have three peaks 
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FIG. 7: The energy E: (top) as function of the fiUing n for 
T = 0.8 and U = 8, 10 and 20, the ED data (4 x 4) are for T = 
and from Ref. l2Rt (middle) as function of the temperature 
T for (7 = 12 and n = 0.5, 0.75 and 1 [(7 = 24 and J = 1/6 
(t-J mode) data are also reported], the qMC data (6x6 and 
8x8) are from Ref. l33 (bottom) as function of the Coulomb 
repulsion [/ for n = 1 and T = 0.8, the ED data (4 x 4) are 
for T = and from Refjpl, the qMCI data (8 x 8) are for 
T = 0.1 and from Ref.jsathe qMC2 data (12 x 12) are for 
T = 1/6 and from Ref.lH 




1.0 




■ 1 ■ 1 ■ 1 


i 1 i 1 i 


1 ' 1 ' 1 ' 1 ' 


0.9 








n = 1 










U = 


0.8 








U = 4 - 


0.7 




I l\ 




U = 8 




(\ 
l\ 




U = 16 " 


0.6 




A 
A 




U = 24 - 






■ A 




U = 32 


0.6 




\ 






0.4 


ii 

i 








0.3 


Ti". 










■ i 








0.2 


ri 










;i 








0.1 


'•i 








0.0 









1 ■ 



9 ^ 10 



FIG. 8: The specific heat C: as function of the temperature 
T for 71 = 0.75 (top) [n = 1 (bottom)] and [/ = 0, 4, 8, 16, 24 
and 32. 



at temperatures of the order half (0.3 0.6) the scales of 
energies Ju, t and U. This is in agreement with the be- 
havior of a pure two- or three- gap system in the canonical 
ensemble. The effect of increasing the Coulomb repulsion 
[/ is a better resolution of the three peaks. In fact, the 
first one moves towards lower temperatures, the second 
one is stable and the third moves to higher temperatures; 
this is in perfect agreement with their origins: l/U, t and 
U, respectively. At half filling nothing changes except 
for the absence of the middle peak: the single occupied 
states do not contribute. It is worth mentioning that the 
kinetic peak appears as the energy levels corresponding 
to different values of the momentum are here discrete as 
in any finite system. No such a peak is present in infi- 
nite systems where the kinetic energy just spreads over 
a band the energy levels coming from the interactions. 

The specific heat C is a really valuable property to 
quantify how good is the t-J model to describe the low 
energy dynamics of the Hubbard model. We have studied 
C as a function of both 4t^/f7 and J in the Hubbard and 
t-J model, respectively. We need a Coulomb repulsion 
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FIG. 9: The specific heat C as function of the temperature T 
for n = 0.75 and [/ = 32 and J = 0.125. 
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FIG. 10: (top) The specific heat C as function of the temper- 
ature T for n = 0.75 and J = 0.0.5 , 0.125, 1/6, 0.25 and 0.5; 
(bottom) The entropy S as function of the temperature T for 
n = 0.75 and J = 0.0.5 , 0.125, 1/6, 0.25 and 0.5. 



V > 24t {At'^/U < to have the possibihty of a faithful 
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FIG. 11: (top) The specific heat C as a function of the tem- 
perature r for i7 = 12 at n = 1. The (7 = 24 and J = 1/6 
{t-J mode) data are also reported. The FTL data (4 x 4) are 
taken from Ref. IsM (bottom) The specific heat C and the 
entropy S as functions of the temperature T for J = 0.15 at 
n = 1. The FTL data (iV = 20) are from Ref.lH 



mapping at low temperatures; only for such values of 
U we manage to sufhciently resolve the three scales of 
energy in the Hubbard model. For instance, at U — 
32 <^ J = 0.125 the mapping can be absolutely trusted 
as shown in Fig. where we report the specific heat for 
n = 0.75. As regards the exchange and kinetic peaks 
the agreement is perfect; obviously the t-J model cannot 
reproduce the Coulomb peak, which is extraneous to its 
dynamics. The absence of the kinetic peak at half filling 
is also perfectly reproduced (see Fig. 1111) . 

In Fig. 1101 (top panel) we report the behavior, in the 
t-J model, of the specific heat as function of the tem- 
perature T for different values of the exchange constant 
J. It is worth mentioning that in the t-J model the 
two peaks in the specific heat (i.e., the exchange and ki- 
netic peaks) can be easily studied separately as they just 
come from the corresponding terms of the Hamiltonian. 
This simple analysis cannot be performed in the Hub- 
bard model where the exchange peak gets contribution 
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from both terms of the Hamiltonian. 

In Fig. ^1 (bottom panel) we report the entropy S of t- 
J model computed by means of the usual thermodynamic 
relations from the specific heat C 



SiT) 



df 



(5.4) 



We can easily put in correspondence the peaks of the 
specific heat [see Fig. ^| (middle panel)] and the change 
in the slope of the entropy. Again, a lowering of the 
exchange energy J helps resolving the scales of energy 
and makes much more visible the difference in the slopes. 

In Fig. ^2 it is reported the specific heat as a function 
of the temperature for the Hubbard model (top panel) 
and the t-J model (bottom panel). For this latter, it 
has been also reported the entropy. The Finite Temper- 
ature Lanczos (FTL) data for the Hubbard model (4 x 4) 
and for the t-J model {N — 20) are taken from Ref. |33 
and Ref. respectively. The agreement is absolutely 
noteworthy on the whole temperature range. Again, the 
discrepancy at low temperatures in the Hubbard model 
case is a consequence of the difference regarding the value 
of the exchange energy J. According to this, we have also 
reported the 2-site results for U = 2A and J — 1/6 {t-J 
model) and obtained the expected agreement. In the t-J 
model case, we have used a value of J twice smaller than 
the one used within the numerical analysis according to 
the boundary conditions we applied to the 2-site system. 
These results show, once more, that the correct and nec- 
essary scales of energies are already present in the 2-site 
system. 



2. The temperature dependence of 
the chemical potential fj, and the double occupancy D 

The temperature dependence of the chemical poten- 
tial /i is given in Fig. E| (top panel), /i has only one 
maximum for U 8 and two maxima and a minimum 
for higher values of U. In the t-J model, the chemi- 
cal potential fi shows a behavior at low temperatures in 
agreement with the one found for the Hubbard model 
(i.e., lim[/^f Ju = J). On the contrary, for high temper- 
atures the divergence is upward instead of downward [see 
Fig. [121 (bottom panel)] . 

The temperature dependence of the chemical potential 
can be explained as follows. In the Hubbard model and 
at zero temperature, there exists a critical value of the 
Coulomb interaction, function of the filling, above which 
the potential energy becomes larger, in absolute value, 
than the kinetic one. This critical value is a decreasing 
function of the filling: lower is the doping higher the dou- 
ble occupancy, bigger the potential energy and smaller, 
in absolute value, the kinetic energy. Then, by increas- 
ing the filling we increase the energy and the chemical 
potential increases. At higher and higher temperatures, 
the system behaves like a free system at zero tempera- 
ture and has a decreasing chemical potential. In partic- 
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FIG. 12: The chemical potential fi: (top) as function of the 
temperature T for n = 0.75 and [/ = 0, 4, 8, 16, 24 and 32; 
(middle) as function of the temperature T for n = 0.75 and 
J = 0.05 , 0.125, 1/6, 0.25 and 0.5; (bottom) as function of 
the temperature T for [/ = 12 and n = 0^ 0.6, 0.7, 0.8 and 
0.9. The FTL data (4 x 4) are from Ref.|H 



ular, for a diverging temperature we have a negatively 
diverging chemical potential (i.e., /i = dE /dn — TdS/dn) 
as the entropy is an increasing function of the filling at 
very high temperatures (i.e., in the almost free case) and 
n < 1. For intermediate temperatures, the chemical po- 
tential has maxima and minima in coincidence with the 
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peaks in the specific heat. Any peak in the specific heat 
marks the temperature at which the system gets the free- 
dom to occupy a new state. The spin (i.e., T ~ Jy) and 
the charge (i.e., T ^ U) peaks reflect in two maxima as 
they signal the availability of triplet states and of dou- 
bly occupied states, respectively. Both states has zero 
kinetic energy. On the contrary, the kinetic (i.e., T ^ t) 
peak signals the availability of further singly occupied 
states which maximize the absolute value of the kinetic 
energy. This reasoning strongly relies on the region of 
filling the figure refers to (i.e., n « 0.75). For other dop- 
ing regions the situation can be quite different owing to 
the behavior in filling of the double occupancy, the en- 
tropy and the positions and existence of the peaks in the 
specific heat. These considerations can also explain the 
quite different behavior we have in the t-J model. At 
n = 0.75 the t-J model is approaching the full-filled sys- 
tem (i.e., n = 1), on the contrary the Hubbard model is 
approaching the half-filled system. At zero temperature 
for n > 0.5, increasing the filling we lowers, in absolute 
value, the kinetic energy as the single electron states, 
which are the only ones with a finite kinetic energy, are 
replaced by the singlet and triplet states: the chemical 
potential is positive. Anyway, the exchange energy can 
effectively lower its value on increasing J. At very high 
temperatures (i.e., in the almost free case), the entropy is 
now a decreasing function of the filling and leads to the 
positive divergence. At intermediate temperatures, the 
reasoning is identical to the one given for the Hubbard 
model except for the obvious absence of the charge peak 
in the t-J model. 

In Fig. (bottom), it is reported the chemical poten- 
tial as function of the temperature T for U — 12 and 
n = 0.5, OA 0.7, 0.8 and 0.9. The FTL data (4 x 4) are 
from Ref . [23. The agreement is really excellent except at 
low temperatures and for the higher values of the filling 
(n = 0.8 and n = 0.9). At low temperatures, owing to 
the finite level spacing the chemical potential has a step 
behavior and, in particular, between n = 0.5 and n = 0.9 
has practically no variation (cfr. FigO. At the higher 
values of the filling (n — 0.8 and n — 0.9), the number of 
2-site states excited at high temperatures is larger than 
those needed to mime the cluster. 

The temperature dependence of the double occupancy 
D at half filling shows a minimum at around T « 1 be- 
fore saturating at the non- interacting value n^/4 for very 
high temperatures^S*^. This minimum coincides with 
the kinetic peak in the specific heat. As already dis- 
cussed for the chemical potential, this peak marks the 
freedom for the system to occupy further singly occupied 
states which have, obviously, zero double occupancy and, 
therefore, lower its total value. 

I 



'2(r, n=l) = ^ 48i2e2/5* (1 _ Se^/^* + e' 

^ ' 8T4(l-f e2/3*)' L ^ 



3. The crossing points 



In Fig. ll3l we report the temperature dependence of the 

specific heat for [/ < 5t at half filling (top panel) and n — 
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FIG. 13: The specific heat C as function of the temperature 
T for n = 1 (left panel) [n = 0.75 (right panel)] and (7 = 0, 
1, 2, 3, 4 and 5. 



0.75 (bottom panel). For these low values of the Coulomb 
repulsion we identify three crossing points/regions where 
the specific heat is almost independent from the value of 
the Coulomb repulsion [/^^L^SiSS. To compute the posi- 
tions and heights of the crossing points we have expanded 
the specific heat C as function of the Coulomb repulsion 
U. We have got 

C(r, U, n) = Co(T, n) + C2{T, n)U^ + O {V^) (5.5a) 
Co(T, n = 1) = 2/3^^2 [1 _ tanh^ (/?<)] (5.5b) 



) 2tT (1 + 20e2^* - 20e^^' - e^''*) -t- 2T^ (l e'^^'f 
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(5.6) 



The temperatures T*(n) at which the crossing point 
can be observed are determined by the equation 
C2{T,n) — 0. At half fiUing and at n = 0.75, this equa- 
tion gives the foUowing results: 



T*in 
T;{n 
T;{n 



T*in 
T*{n 



1) = 0.499 
1) 0.997 
1) = 2.058 



0.75) 0.951 
0.75) 2.024 



Co =^ 0.563 
Co =s 0.841 
Co = 0.376 



Co = 0.523 
Co = 0.342 



(5.7a) 
(5.7b) 
(5.7c) 



(5.8a) 
(5.8b) 



VoUhardt found an approximate formula for Co and T*, 
as function of the dimensionality of the system, for the 
higher temperature crossing point at half filling"^^. The 
values that this formula gives for d = 1, the dimension- 
ality of the 2-site system within periodic boundary con- 
ditions, are close, as regards Co, to the exact ones of 
Eq. f57c|) . 

It is worth noting that higher is the temperature nar- 
rower is the crossing region; ranging from quite wide re- 
gions for the two lower temperatures to a very sharp point 
for the higher temperature. Also the chemical potential 
(Fig.El(top panel) at n « 0.275, 0.75, 1.25, 1.725) and the 
double occupancy show quite clear crossing points in (no 
dependence on) temperature. 
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E. Charge, spin and pair correlations 



FIG. 14: The charge x? (top) and spin (bottom) correla- 
tion functions as functions of 17 at n = 1 and for T = 0, 0.1, 
0.5 and 1. 



The charge (i.e., Xc)i spi'^ (i-^-j x" ) and pair (i.e., d) 
correlation functions contain information regarding the 
spatial distributions of the corresponding quantities. Ob- 
viously, we can only consider the quantum fluctuations 
in the paramagnetic state, which is the only admissible 
equilibrium state on a finite cluster. 

The charge correlation function x" is shown in Fig. 1141 
(top) as function of the interaction U for different values 
of the temperature T and n — \. In the non- interacting 
case {[] — 0) and at zero temperature we have 



if n < 1 
if n > 1 



(5.9) 



as only the singlet states with one electron {n <V) [one 
hole (n > 1)] per site contribute. In fact, these are the 
states which lower the most the internal energy as they 
maximize the absolute value of the kinetic one. In the 
strongly interacting limit {JJ — > oo) and at zero temper- 




2 4 6 



FIG. 15: Comparison between Jy and its limiting value At^ /U 
as functions of U . 
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FIG. 16: The pair correlation function d as function of T at 
n = 1 and for [/ = 0, 2, 4, 8, 12 and 16. 



ature we have 

{-n2 if n < 0.5 

if 0.5<n< 1.5 (5.10) 
-{2~nf if n> 1.5 

as only the states with a single electron {n < 0.5) [a 
single hole (n > 1.5)] or only with one electron per site 
(0.5 < n < 1) [one hole per site (1 < n < 1.5)] contribute. 
In fact, no double occupancy is allowed and there is no 
gain in the kinetic energy if the singlet states are used as 
J(7 0. At intermediate values of the coupling U the 
double occupied states play a relevant role and we found 
results between the limiting cases reported above: on in- 
creasing the coupling U the double occupancy diminishes 
and consequently the states with one electron per site in- 
crease their contribution and the value of x" (see Fig. 
(top)). On increasing the temperature T, more and more 
states become available and the correlation tends to its 
reducible part, n^. 

The spin correlation function is reported in Fig. 1141 
(bottom) for the same parameters chosen for In the 
non-interacting case {U — 0), we have Xs = Xc- i^i 
absence of interaction the charge and spin behave co- 
herently. On increasing the interaction potential U the 
singlet states are favored and the spin correlations get 
enhanced {Jjj gets closer and closer to At^ /U which is 
the scale of energy of the spin excitations, see Fig. [T^ . 
Then, further increasing U the value of Ju tends to zero 
and the spin correlations get suppressed except at zero 
temperature [see Fig. 1141 fbottoml] . A rapid suppression 
of the spin fluctuations can be caused also by the incre- 
ment of the temperature. 

The pair correlation function d is reported in Fig. llBl as 
function of the temperature T at n = 1 and for ?7 = 0, 2, 
4, 8, 12 and 16. Obviously, d decreases on increasing the 
Coulomb interaction U and is maximum at half-filling 
where the singlet state is the favorite one. On increasing 
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FIG. 17: The charge Xc(0,0) (top) and spin (0, 0) (bottom) 
susceptibilities as functions of T at n = 1 for U = 0, 2, 4, 8, 
12 and 16. 



the temperature T (see Fig. ll6|l . more and more states be- 
come available and the pair correlation function d tends 
to zero. 



F. Charge and spin susceptibilities 



by 



The charge and spin susceptibilities Xfiik,uj) are given 



(5.11) 



where no summation is implied by repeated indices. By 
means of the expression of the causal Green's function 
given in Sec. lIVB 1 al ls immediate to compute the follow- 
ing static susceptibilities through the spectral theorem 



iXc(fc,0) = -C,",-^d 



ryot 



1 U 



(5.12) 
(5.13) 



where d and x" are the first-neighbor pair and spin cor- 
relation functions, respectively. 

In Fig. El the charge Xc{k,0) and spin Xs{k,0) sus- 
ceptibilities are reported as function of the temperature 
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FIG. 18: The compressibility lunction of T at n = 1 



T for different values of the Coulomb repulsion IJ and 
71 = 1. As it results clear by Eqs. 15.121 the two sus- 
ceptibilities result identical in the non-interacting case 
[/ = 0: once more charge and spin behave coherently in 
absence of interaction. On increasing the Coulomb repul- 
sion, the two susceptibilities behave in opposite manners: 
the charge susceptibility is strongly suppressed, in par- 
ticular at half-filling where a charge gap develops in the 
strongly interacting limit {U — *■ cx)); the spin susceptibil- 
ity is greatly enhanced, in particular at half-filling where 
the singlet state is the favorite one. The behavior of 
the two susceptibilities as functions of the temperature 
deserve special attention as it is directly connected with 
the formation of a gap in the corresponding channel. The 
spin susceptibility shows the typical paramagnetic behav- 
ior: plateau at low temperatures and Curie tail at high 
temperatures. The charge susceptibility instead shows 
the presence of a gap in the strongly interacting limit 
([/ oo). The small upturn at low temperatures and 
medium-high values of the Coulomb repulsion V is due 
to the finite level-spacing of the system. In particular, it 
arises from the intra-subband gaps characteristic of the 
scale of energy Jj/. 

In the t-J model, we simply have 



iXc(fc,0) = -C7fi 
txs (fc,0) = -Cfi 



1 t 
3 



(5.14) 
(5.15) 



G. Thermal compressibility 

The thermal compressibility is defined as 
1 dn 



(5.16) 



By using some general quantum statistical relations, 
which can be established between the particle density 



n and the chemical potential /j^, we can express the 
thermal compressibility in terms of the density-density 
correlation function 



1 



2L> 



2 



) + (x? - -^)] 



(5.17) 



In Fig. 1181 the thermal compressibility k is reported 
as function of the temperature T for different values of 
the Coulomb repulsion V and n — \. The system is 
completely incompressible k = when no more particles 
are allowed to enter the system: the chemical potential, 
which is a measure of the energy necessary to insert a new 
particle in the system, diverges. Obviously, the system 
is extremely eager to accept particles at very low fillings, 
in order to increase in absolute value the kinetic energy 
and lower the total one, and absolutely incompressible at 
n = 2, when all the quantum states are filled. At zero 
temperature, the Coulomb repulsion makes incompress- 
ible also the states with commensurate fractional fillings 
n = 0.5 and n — 1.5. At half-filling, on increasing the 
Coulomb repulsion, the compressibility is rapidly sup- 
pressed according to the very high price in energy that 
should be paid to add one particle to the singlet state 
which is the one favored by the Coulomb repulsion. The 
most relevant feature for this property is the presence of a 
well-defined peak when it is plotted versus temperature. 
A finite, but low in comparison with the actual value of 
the Coulomb repulsion, temperature permits to overcome 
the suppression related to the formation of a charge gap. 
A further increment of the temperature makes available 
more and more states and the system is driven back to 
be incompressible. 

It is possible to write the compressibility for the Hub- 
bard and t-J models in an identical way 



Tn- 



(File 



(5.18) 



where Tuc is the zero frequency function in the charge 
channel and its ergodic value is exactly n^. According 
to this, the compressibility is a direct measure of the 
ergodicity of the charge dynamics. 



H. Optical conductivity 

In the framework of the linear response and by us- 
ing the Ward-Takahashi identities'^ , which relate the 
current-current propagator to the charge-charge one by 
exploiting the charge conservation, the optical conduc- 
tivity is given by 



ai{uj) = Du]5{uj) 
with the Drude weight Dw given by 



(5.19) 



(5.20) 



For the 2-site system, the incoherent part of optical con- 
ductivity is zero as no contribution comes from the imag- 
inary part of retarded current-current propagator. 
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FIG. 19: The Drude weight as a function of n for T — 1 



In Fig. m the Drude weight D^, normahzed by 
27re^a^, is reported as function of the fiUing n for dif- 
ferent values of the Coulomb repulsion U and T = 1. 
At half-filling, the Coulomb repulsion tends to suppress 
the Drude weight and drives the system to be insulating. 
However, the Drude weight vanishes only in the limit U 
infinite. Higher the temperature more states with no con- 
tribution to the kinetic energy result available and lower 
is the Drude weight. 
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VI. ERGODICITY 

One of the main issues on which this manuscript wishes 
to draw attention is the non-ergodicity of the charge and 
spin dynamics in the two-site Hubbard model. The non- 
ergodic dynamics in this system is due to the finite num- 
ber of degenerate states available to the system when 
dealing with finite temperatures and/or incommensurate 
fillings and/or no interaction (cfr. Tableland Ea. lA12| . 
The finite level spacing confines the system to degen- 
erate or non-degenerate states, according to the filling 
and interaction strength, making the dynamic ergodic 
or non-ergodic according to the degree of degeneracy of 
the ground state. The temperature then opens up the 
possibility for quite complicate mixtures, only possible 
according to our choice to work in the grand-canonical 
ensemble in order to get results comparable with larger 
systems for which the level spacing is thinner. The cou- 
pling to an heat and a particle reservoirs have given us the 
possibility to simulate a continuous tuning of the filling 
and temperature on a finite system (with a finite degree 
of freedom) as it is possible only in a bulk system with 
infinite degrees of freedom. 

The zero-frequency constants for charge Fuc and spin 
Fiis channels have been plotted as functions of the fill- 
ing n for different values of the Coulomb repulsion U and 
T = 0.1 and as functions of the interaction U for differ- 
ent values of temperature T and n = 1 in Figs. [201 and 
im Their ergodic values are and 0, respectively. Fuc 



FIG. 20: The charge zero frequency constant Fuc diminished 
of n^: (top) as function of n for T = 0.1 and (7 = 0, 2, 4, 8, 
12 and 16; (bottom) as function of f/ at n = 1 and for T = 0, 
0.1, 0.5 and 1. 

assumes its ergodic value only at zero temperature and 
commensurate fillings; in particular, only at integer com- 
mensurate fillings for any interaction strength and also 
at fractional ones for a finite interaction strength. On 
the other hand, Fn^, assumes its ergodic value only at 
zero temperature and integer commensurate fillings for 
any interaction strength. At half-filling, while the charge 
dynamics is constrained to be ergodic in the strongly in- 
teracting limit (U oo) by the formation of a charge 
gap, the spin dynamics gets more and more non-ergodic 
as larger the interaction strength is. The states that are 
left available in this condition have a quite different be- 
havior once probed with respect to the charge or spin 
response. 

VII. CONCLUSIONS 

We have studied the 2-site Hubbard and t-J models by 
means of the Green's function and equations of motion 
formalism. The main results can be so summarized: 

• We have got a complete basis of eigenoperators for 
the fermionic and bosonic sectors which could be 
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FIG. 21: The spin zero frequency constant Fus'. (top) as 
function of n for T = 0.1 and U = 0, 2, 4, 8, 12 and 16; 
(bottom) as function of f/ at n = 1 and for T = 0, 0.1, 0.5 
and 1. 



report that the exact results of the 2-site system 
managed to reproduce the Mott-Hubbard MIT of 
the bulk, the behavior of some local quantities (i.e., 
the double occupancy, the local magnetic moment 
and the kinetic energy) and of some thermody- 
namic quantities (i.e., the energy, the specific heat 
and the entropy) of larger clusters. This shows that 
the necessary energy scales are already present in 
the 2-site system. According to this, as already said 
in the first point above, we strongly believe that the 
operatorial basis that exactly solves this system can 
also give excellent results if used for the bulk. 

• The study of the specific heat has given many valu- 
able information regarding: the scales of energy 
present in the system, their origin, interaction and 
possibility of resolution, the range of parameters 
for which the t-J model faithfully reproduces the 
low energy /temperature behavior of the Hubbard 
model, the explanation for the temperature depen- 
dence of the local properties, the existence of cross- 
ing points. 

• We have shown how relevant is the determination 
of the zero-frequency constants in order to correctly 
compute the bosonic Green's functions. In particu- 
lar, we have shown that, for the 2-site system, they 
assume values very far from the ergodic ones. 
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used to get a controlled approximation in the study 
of the lattice case by means of any approximation 
that strongly relies on the choice of the basic field. 

• We have identified the eigenoperator responsible for 
the appearance of the exchange scale of energy J. 

• We have illustrated, once moreil, as the local alge- 
bra constrains can properly fix the representation 
and easily give the values of the bosonic correlations 
that appear as internal parameters in the fermionic 
sector and of the zero frequency functions that ap- 
pear in the bosonic sector. This also permits, while 
studying the fermionic sector, to avoid the opening 
of the bosonic one (i.e., the charge, the spin and 
pair channels) and all the heavy calculations re- 
quired to solve it. 

• We have explored the possibility of the 2-site sys- 
tems to mime the behavior of larger clusters, as 
regards some of their physical properties, by using 
an higher temperature and by exploiting the qual- 
itative properties of the level spacing. We have to 



APPENDIX A: THE EIGENPROBLEM 

The eigenproblem of a given grand canonical Hamil- 
tonian H is solved once the latter is diagonalized on the 
Fock space of the system under study, that is 

H \n) = En (Ala) 
N \n) = Nn \n) (Alb) 

where N = n(i) is the total number operator and \n) 
is a complete orthonormal basis. 

Once the eigenproblem is solved, we can compute the 
thermal average of an operator <& by means of the follow- 
ing expression 

n 

where Z — J2n e^'^-^" is the grand canonical partition 
function and (3 is the inverse temperature. 
If a finite minimal energy i^min exists then 

^z'^'^^^j^ E (A3) 

n\E„=E^-nin 
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where M is the number of eigenstates \n) with En = 
E ■ 

On the basis of the knowledge of the set of eigenstates 
and eigenvalues of H l|Ala|) it is possible to derive the ex- 
pressions for Green's functions and correlation functions. 
Let be a field operator in the Hciscnberg scheme 



where 

^^'^(iJ) = (n|V'(i)|m)(TO|V't(j)|n) 



(All) 



and the zero frequency constant F^^t (i, j) has the follow- 
ing representation 



i'^*-0(i)e we do not specify the ^ (■ •\ — }_ \^ -/3S„ ^n.m/. 
r bosonic. nUhii^ (ih(i) ca,n be. for in- V"AtU,Jj - 7 <^ ^v^tliJJ 



nature, fermionic or bosonic, ofip^i) {ip{i) can be, for in 
stance, either c{i) or n(i)). By considering the two-time 
thermodynamic Green's functions^i*^'" , let us define the 
causal function 



(A12) 



-rje{t,~U){i:^{j)^{i)) (A4) 
the retarded and advanced functions 

G^■^,^'^\^,J) = ±6 [± iU - t,)]{[m,^Hj)],) (A5) 
and the correlation function 

C^i.f{^,J) = {m^\J)) (A6) 

Here 77 = ±1; usually, it is convenient to take 77 = 1 
(77 = —1) for a fermionic (bosonic) field ^p (i-e., for 
a composite field constituted of an odd (even) num- 
ber of original fields) in order to exploit the canonical 
(anti)commutation relations of the constituting original 
fields; but, in principle, both choices are possible. Ac- 
cordingly, we define 

_j{AB} = AB + BA for ,7 = 1 
" \ [A, B]=AB- BA for 77 = -1 



By using the definition of thermal average ljA2|l it is pos- 
sible to derive the following expressions for the Green's 
functions (|A4I IA5|) and correlation functions l|A6|) in 
terms of eigenstates and eigenvalues 



G^^'!-^(i,j,w) = r^^t(i,j) 



1 



{I + f])V - - iTT{l - T^)5{UJ) 



- y 

z ^ 

E„^E„ 



A 



,A;-(i,j)e-^^" 



UJ + En - Em + iS 



En — Em — i<5 



(A8) 



Gj;,r)(i,j,-) = r^^t 



(U) 



1 



V 



- y 

z 



AS(i,j)(e-^^" + 



77e 



n , m 



En — E„i ± id 



(A9) 



E„=E„ 



It is worth noticing that, in this presentation (i.e., in 
terms of eigenstates and eigenvalues), the Green's func- 
tions and correlation functions are fully determined up 
to the value of the chemical potential, which is present 
in the expressions of the eigenvalues. As usual, the com- 
putation of the chemical potential requires the inversion 
of the expression, in terms of eigenstates and eigenvalues, 
for the number of particle per site (|A42() . In the main text 
instead, we presented expressions for the same quantities 
(Green's and correlation functions) in terms of eigenen- 
ergies of eigenoperators and of correlation functions of 
these latter. In this case, more parameters appeared 
(A and p in the fermionic sector and x"i d, 7 and the 
zero frequency constant Fn^ in the bosonic sector) that 
have been computed self-consistently as we usually do for 
the chemical potential. The two presentations, although 
equivalent (i.e., they obviously lead to the same results), 
require two different self-consistent procedures to come to 
the computation of the physical quantities. The reason 
of this occurrence resides in the level of knowledge of the 
representation we have in the two cases. In the first case, 
we have the full knowledge of the states of the system and 
we should only fix their occupancy with respect to the 
average number of particle per site we wish to fix; in the 
second case, as we have no direct knowledge of the states 
before fixing the counting we have to reduce the Hilbert 
space to the correct one (i.e., the one of the system under 
analysis), that is, we have to impose constraints in order 
to select only those states that enjoy the correct symme- 
try properties (e.g., we have to discard states with site 
double occupied by electrons with parallel spins) . At any 
rate, we want to emphasize once more that all the results 
obtained in this paper by means of the Green's function 
formalism exactly coincide, as it should be, with those 
computable by means of the thermal averages. 

The eigenstates of the systems under study are given 
by linear combinations of the vectors spanning their Pock 
spaces. These vectors can be displayed as \a, h) where a 
and h denote the occupancy of each site, respectively. In 
particular, we have for an empty site, "f (|) for a single 
occupied site by a spin- up (spin-down) electron and "fj, 
for a double occupied site. This latter is not allowed in 
the t-J model (see Tab.Pl. 



C^V'i"(i'j'^) = 27rF^^t(i,j)^(t^) 



27r 



+ ^ E e-p''-A';;,(\rm^ + En-Em) (aio) 



Er,^E„ 



The Hubbard model 



The eigenstates and eigenenergies of the 2-site Hub- 
bard model are reported in Tab. ^ The coefficients 
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7i[|T,o)-|o,T)] 


-^i + 2t 
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1 
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3 


75[|i,0)-|0,i)] 


-fi + 2t 


1 


1 
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;^[|T,0) + |0,T)] 


-fi - 2t 


1 


1 

2 


5 


;^[U,o) + |o,i)] 


-fi - 2t 


1 


_1 
2 


6 


IT,T) 




2 


1 


7 




-2^i 


2 


-1 


8 


;^[|T,i) + li,T)] 


-2^i 


2 





9 


ai [iT,i)-ii,T)]-a2[in,o) + |o,n)i 


-2^4 - 4:Ju 


2 





10 


75[|Ti,0)-|0,n)] 


-2fi + U 


2 





11 


"2 [iT,i) -u,T)] + ai [in,o) + |o,n)] 


-2p + f/ + 4Jc7 


2 





12 


73[in,T) + iT,n>] 


-3At + 2t + U 
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1 
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13 


;^[in,i) + ii,n)] 


-3^ + 2t + i7 
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1 
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14 


75[|Ti,T)- IT,n)] 


-3At - 2t + t/ 
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1 
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15 


75[in,i)-ii,n)] 


-3^l -2t + U 


3 


1 
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16 


in,n) 


+ 2U 


4 






TABLE I: The eigenstates and eigenenergies of the 2-site Hubbard and t-J (first nine with qi = l/\/2, Q2 = and Jjj = J) 
models. In the last column is reported the corresponding eigenvalue of Sz that, together with E„ the eigenvalue of the 
Hamiltonian, completely characterized the states (i.e., the presence of an external magnetic field will completely lift the 
degeneracy) . 



ai and 02 are determined by the orthonormality of the where 
eigenstates and have the following expressions: 

a, = (U + ^Ju)^ (A13) S;," = I (ai,2 + A4^6 + ^As.s + a?A3,9 



2ty^ + 2 ( '^zAs.ll + OilAg^i4 + -Aio,12 + tti Aii_i4 



a2 = = (A14) 

\/{U + AJuf + 16^2 



(A19) 



(A20) 



We also have a\+ a\ — 1/2. 

Expressions (|X8)l . (|X9)) and l|A10(l show that the 

Green's functions and the correlation functions are com- ^n,m 1 ^ 1^ 2^ 

pletely determined (up to the value of the chemical po- 5^5]; ~ 2 I i-'' ^ ^ 2 ^'^ ^ ^'^ 

tential) once the matrices A^'^(i, j) are known. We here ^ / ^ \ 

give the results for some relevant operators (i, j = 1, 2; ~'~ 2 V '^^As.ii + a2A9,i2 + 2 Aio,i4 + Aii^i2 j 

Operator ^(i) 

A"'" (1, 1) - B"'" + C"'™ (A15) 

2) = + C^"^'™ (A16) B";^ = - (^Ai,3 + A5,7 + - A4^8 + a? A2,9 

^ ^ ^ +- U2^2,ii+a2A9,i5 + -Aio,i3 + a2An,i5 
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C^"'™ = i (^Ai,5 + A3,7 + ^A2,8 + af A4,! 

+ i f aiA4,ii + aiAg^ia + ^Aio.is + cti An, 



13 

(A22) 



Operator r]{i) 






A"'™ (1,1) = 




(A23) 


A"'" (1,2) = 


-B t + ^ ■ t 


(A24) 


^"'"1(1,1) = 


B t + ^ t 


(A25) 


yl"''"(l,2) = 

VlVi 


-B t + ^ t 


(A26) 


where 







C'^iT = 4 (A2,4 + A3,5 + A4,2 + A5,3) 

+ ^ (Ai2,l4 + Ai3,i5 + Ai4,l2 + Ai5,i3) (A34) 



D2;r = 2 (^2^9,10 + "2^10,9 + a?Aio,ii + a?An,io) 



Operator 713(1) 



n.7n 
"ana 



4n,m /-I r>n,m 



where 



n3»3 nans 
n3"3 n3"3 



(A35) 



(A36) 
(A37) 



B"'™ = i ( A7,i5 + Ai5,i6 + ^A5,io + A3,ii 



+ i (alAa^g + Q!^Aii,i2 + ^Ag.M + afAg, 



C"'"; = J ( A7,13 + Ai3,16 + ^A3,10 + a?A5,il 
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(A27) 



+ ^ ( ajAs.g + q;2Aii,i4 + :r A8,i2 + a^Ag^^ 



(A28) 



B^;™ = - (A2,2 + A3,3 + A4.4 + A5,5 + 4A6,6) 

+ i (4A7,7 + Ai2,12 + Ai3,i3 + Ai4,14 + Aig^is) 



r\n,7n 
nsna 



(A38) 
(A39) 



2 (alAg,!! + aiAii,8 + a^Ag.g + afAg^g) 



(A40) 



^n,m ^ 1 / ^ ^ ^ 1 ^ ^ a? A2,ii 

2 y 2 



+ ^ f aiA2,9 + Q!2Aii,i3 + ^A8,15 + Q!lAg,i3 



(A29) 



C"'™ = ^ ( A6,12 + Ai2,16 + A2,10 + af A4,ll 



i faiA4,9 + a2Aii,i5 + ^A8,i3 + a? Ag^ig 



Operator n(i) = c^{i)c{i) 

^nn l-"-, -"-J — -°nn + '-'nn + ^nn 
^nn (1, -^) = -Dnn " (^nn " -^nn 

where 

= J (A2,2 + A3,3 + A4,4 + A5,5 + 4A6,6) 
+ (A7,7 + A8,8 + Ag^g + A1040 + An^n) 



(A30) 



(A31) 
(A32) 



2. The t-J model 

The eigenstates and eigenenergies of the 2-site t-J 
model coincides with the first nine of the Hubbard model 
with ai = I/V2, a2 ~ and Ju = J- Actually, in the 
strong coupling regime (i.e., U ^ t) we have ai — > 1/V2, 
a2 — > and Ju — * At^ /U. In particular, this latter value 
{At^ /U) is the one wc get in the strong coupling regime 
when deriving the t-J model from the Hubbard one. All 
other states of the Hubbard model can not be realized 
in the t-J model owing to the exclusion of the double 
occupancies. 

Following the same prescription (only first nine states, 
ai = l/-\/2, a2 = and Jjj = J), it is possible to get 
A^'™(i,j), A»ir(i,j) and (i,j) form the correspond- 
ing expressions given for the Hubbard model (yl^^™(i, j), 
^nr(i.j) and A^J^™ (i,j), respectively). 



3. The properties 



+ ^ (9Ai2,i2 + 9Ai3,i3 + 9Ai4,i4 + 9Ai5,i5 + 16Ai6,i6) 



(A33) 



Given the expressions of A^';^"(i,j), we can now com- 
pute the relevant physical quantities. 
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a. The Hubbard model 

The partition function is given by 

Z = l + 2e-'5^= + 2e-'^^* + Se"''^'^ + e''^^^ + e-^^^' 



Spin (Xs ), charge (x") and pair (d) correlation func- 
tions are given by 

X« = I (e-'^^^ - 2a2e-/5s^ - 2a2e-/5S" ) (A48) 



+ e 



(A41) 



The chemical potential \x can be computed inverting 
the following expression for the particle number per site 



= ^ (Sc-/^^^ + 2af e-'^^^ + 2a2e-'5^^ 

Zj 

+4e-/3Si2+4e-/3iJi4+4e-/3iJia) 



(A49) 



+e-/^^" + Se-'^^- + 3e-^^" + 2e-^^-) (A42) ^ ^ J_ (2«2^-/3i^. _ ^-0e.. + 2a?e-''^") (A50) 



This expression determines the chemical potential as a 
function of n, T and U . 

The self-consistent parameter A and p are given by 

A = ^ (e-'^^^ - e-'^^* - e-^^^^ + e'^^^^^) (A43) 



^2e-/3Si4 + 2e-^^''>) (A44) 
The correlation functions Cf2 and C"^ are given by 

C- = ^ f e-'^^^ - e-'^^" ) (A45) 



^12 - ^ (e 



+e- 



-/3Bi2 _ p-/3Bi 



(A46) 



The double occupancy per site is given by 

2Z ^ 



+2e-^^i* + 2e-'^^i^) 



(A47) 



2Z 

The zero-frequency constant Fno = rnn(l,2) and 
Tiis = r„3„3(l, 2) are given by 

Tno = — (c-^^^ + c-^^^ + 6c-'3^« + 2c-^^^» + 2e-'5^^° 



+2e 



9e 



-/3Bi6 



(A51) 



2Z 



(A52) 



4. The t-J model 



The partition function is given by 

Z=l + 2e-'^^' + 2e-'^^* + Se"''^'' + e'^'^' (A53) 

The chemical potential fj, can be computed through the 
following expression for the particle number per site 

n=he-^^'+ e-^^- + Se''^^^ + e"^^^ ) ( A54) 
z 

which can be inverted and give 



(2n - 1) cosh(2/3i) + J (2n - 1)2 cosh^(2l3t) + n(l - n)(3 + e^^^) 
^ = ^^^ (l-n)(3 + eW) 



I 

The self-consistent parameter C^i and x" are given by Spin (x") and charge (x" ) correlation functions are 

given by 

C?i = ^(e-'^^^-e-''^^) (A56) = ^ (e-^E, _ ^-0e,^ (A58) 



x2 = §(3e-^^^-e-/'^«) 



(A57) ;^« = ^(3e-'^^^+e-/'^«) 



(A59) 
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The zero-frequency constant Fno = rnn(l,2) and 
rii3 = r„3„3 (1, 2) are given by 

rno = ^ (e-^-^ + e-^-^ + Ge"'^-^ + 2^-^-^) = ^ (e"^-^ + e"^-^ + 4e-^-«) (A61) 

(A60) 
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